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Abstract 

If more than one curvaton dominate the Universe at different epochs from each 
other, curvature perturbations can be temporarily enhanced to a value much larger 
than the observed one 10~^. The traces of the enhancement may be left as higher 
order correlation functions, that is, as non-Gaussianity, the stochastic gravitational 
waves that are sourced by scalar-scalar mode couplings, as well as the primordial 
black holes that are formed by the gravitational collapse of the enhanced curvature 
perturbations. We first confirm that such a temporal enhancement indeed occurs 
by solving the linearized perturbation equations both numerically and analytically. 
We then derive an analytic expression of the full-order curvature perturbation which 
does not rely on the frequently used sudden decay approximation and is exact on 
super-horizon scales. By using this analytic formula, we provide expressions of the 
non-linearity parameters /nL) tnl and 5nl- If both two curvatons contribute to the 
final curvature perturbations, the strong non-Gaussianity appears in the trispectrum 
rather than in the bispectrum. We also find a unique consistency relation between 
tnl and gfNL without /nl- By using the second-order perturbation theory, we nu- 
merically show that the spectrum of the induced gravitational waves has a plateau 
corresponding to duration of the enhancement and such gravitational waves can be 
probed by ultimate-DECIGO and space-based atomic interferometers. We finally 
calculate the abundance of the primordial black holes and put a constraint on the 
amplitude of the enhanced curvature perturbations. 



1 Introduction 



Cosmic observations [I] are now strongly supporting the idea of primordial inflation pHl]. 
According to the inflationary scenario, the primordial perturbations, which are seeds for 
the inhomogeneous structure of our Universe, are created from the field fluctuations that 
are generated quantum mechanically during inflation [6H9]. In the simplest inflationary 
scenario, the fluctuations of the inflaton, where energy drives accelerated expansion, turn 
into the curvature perturbations when its wavelength becomes larger than the Hubble 
length during inflation. Once generated, the curvature perturbation remains constant and 
this is the primordial perturbation that can be compared with CMB observations. 

This simple scenario, however, may not be the one that actually happened in our 
Universe. For example, in the curvaton scenario [T0HT2] . a curvaton, a scalar field other 
than the inflaton, plays a role of creating the primordial perturbation. In this scenario, 
the curvaton perturbations turn into the curvature perturbation after inflation when the 
curvaton decays into the radiation. There are many other models in which the primordial 
perturbations are created after inflation such as inhomogeneous end of inflation [TSHIZ], 
modulated reheating scenario [18l[T9] etc.. A common feature among all these models is 
that the curvature perturbation, once generated from zero remains constant outside the 
horizon. 

Although these nonstandard scenarios are already contrived, the real Universe may 
have evolved even more complicated manner, say, with multiple curvatons dominating at 
different epochs. Such a scenario has been studied in [20] and [21] , the former focusing on 
the power spectrum and the latter on the bispectrum. In [21], although not emphasized, it 
was implicitly shown that if the two curvatons dominate the Universe at different epochs 
from each other, the curvature fluctuations may evolve in a dramatically different way 
than the standard case. That is, they can grow to an amplitude much larger than the 
observed value, 10~^, when the first curvaton dominates and decays, and then they are 
moderated to the observed amplitude when the second curvaton dominates and decays. 
Thus in this scenario the curvature perturbation can be temporarily enhanced. 

At first glance, this temporal enhancement seems to have little effect on observables 
since it must occur, if at all, much before the big-bang nucleosynthesis, not to mention 
the observation time. But this is not the case. We can provide at least three possible 
interesting consequences from this effect. 

The first possibility is that, as pointed out in [21], large possibly detectable non- 
Gaussian perturbations can be generated. If we denote by Cmax the maximum amplitude 
of the curvature perturbation when it is enhanced, the so-called /nl parameter is given by 
/nl — Cmax/10~^ ^ 1, under the assumption that the curvature perturbation is sourced 
only by the first decaying curvaton fluctuations. For example, if Cmax = 10~^, then we get 
/nl ^ 100. 

The second possibility is generation of stochastic gravitational waves(GWs) whose peak 
frequency can fall into a range of the GW detectors. It is well known that at second order 
in perturbation the scalar-scalar coupling can source the GWs, most efficiently when the 
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scalar mode reenters the horizon [22] - [2^ . Since the amphtude of GWs is proportional to 
the square of the curvature perturbations, we expect that large amplitude of GWs can be 
generated at the horizon scales when the curvature perturbation is being enhanced. 

The third possibility is a formation of primordial black holes (PBHs) [251 - I27] . If Cmax 
is very large, the perturbation mode which reenters the horizon when the curvature per- 
turbation is being enhanced may undergo a gravitational collapse to form a black hole 
^^1 . Since the abundance of PBHs on various masses are tightly constrained from cosmic 
observations (see [28] for the latest results), we can limit a range of the parameters of the 
two curvaton model by using such constraints. 

The present paper aims to take up this temporal enhancement of the curvature per- 
turbation in two curvaton model, to provide a detailed analysis of the generation and 
evolution of fluctuations and to discuss three observational implications mentioned above. 



2 Basic picture 

In this paper, we are interested in a situation where both curvatons dominate the Universe 
at different epochs. Before going to the detailed analysis of the scenario, let us first explain 
four non-trivial assumptions needed for the designing scenario to work, along with the basic 
history of the universe under consideration. 

The first assumption we will take in this paper is that, other than the infiaton, there 
are two light free scalar fields (we call them cxi-field and (T2-field, respectively.) having 
VEVs smaller than the Planck scale in the early universe. By light, we mean that masses 
of both curvatons are much smaller than the expansion rate of the Universe when it is 
reheated by the decay of the infiaton, mi, m2 <^ H^^h with nii and m2 being masses of the 
two curvatons and ifreh the Hubble parameter at the time of reheating. This condition 
means that both curvatons are almost massless during inflation. It is well known that such 
a light scalar field aquires classical fluctuation of order ifinf/27r on super- horizon scales. 
The requirement of VEVs smaller than the Planck scale is to avoid the second inflation 
caused by a curvaton. Although this requirement is not essential, we take it for simplicity. 

After inflation ends, the infiaton oscillates around the minimum of the potential and 
finally decays into radiation to complete reheating. Then the radiation dominated universe 
starts. At this stage, both curvatons are still subdominant. Therefore, those fluctuations 
have not yet contributed to the curvature perturbation. 

The next assumption is that both curvatons start to ocsillate during this radiation 
dominated epoch. If the secondly decaying curvaton starts to oscillate after the first 
decaying curvaton decays, the temporal enhancement of the curvature perturbation, which 
is the main focus of this paper, does not occur. Therefore, we do not consider such a case. 

"^^Some papers, such as [291135] . also considered the generation of large amplitude of the curvature 
perturbations on particular scales and the PBH formation. Contrary to the two curvaton model considered 
in this paper, the curvature perturbations in those models do not show the temporal enhancement on 
super-horizon scales. 
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Since the ocsillating free scalar field can be treated as non-relativistic particles, we 
can assume that the universe consists of the radiation coming from inflaton, two non- 
relativistic particles that are not interacting. During this epoch, the fraction of energy 
densities of both curvatons grows in proportion to the scale factor while the ratio between 
energy densities of the two curvatons stays constant. We call this era epoch A. 

Without a loss of generality, we can assume that ai-field decays first and cr2-field 
decays later, Fi ^ Then, the third assumption is that the energy density of ai- 
field in epoch A is much larger than that of cT2-field and there is a period when ai-field 
dominates the universe. We call such a period epoch B. Since the cxi-field is dominating 
the universe, it is this epoch when the o"i-field perturbation is mostly converted to the 
curvature perturbation. By exactly the same reasoning as the single curvaton case, the 
curvature perturbation at this epoch is given by ~ 5i = Spi/pi (precise definition of 
the curvature perturbation and regorous calculations will be given later.). Since the o"2- 
field is subdominant during this epoch, keeping its energy fraction constant, the (T2-field 
fluctuations have not been converted to the curvature perturbations yet. When the Hubble 
parameter becomes Fi, the cxi-field decays into the radition and the universe is again 
dominated by the radiation. We call a period dominated by such radiation epoch C. 

Our last assumption is that the (T2-field finally dominates the universe before it decays. 
In other words, r2 should be small enough to allow the cr2-field to dominate the universe. 
We call a period dominated by the (T2-field epoch D. During this epoch, since the universe 
evolves like a matter dominated universe, the radiation generated from the decay of the 
(Ji-field is diluted. Because of this dilution, the curvature perturbation coming from the 
radiation perturbation is reduced by a factor fi^ <C 1, where Qr is a fraction of the 
radiation energy density to the total one. Meanwhile, the (T2-field contributes to the 
curvature perturbation by 62- When the Hubble parameter becomes equal to r2, the (T2- 
field decays into radiation and the Universe is again dominated by the radiation. Since 
there are no isocurvature perturbations any more, the curvature perturbation remains 
constant and this should be regarded as the final perturbation that can be compared with 
the cosmological observations such as CMB. Thus, the primordial perturbation at the 
linear order is estimated as 



where fi^ must be evaluated at the time when (T2-field decays. From observations, we know 
that this is about 10~^. 

From those arguments, we find that the curvature perturbation evolves from zero 
to 61 at epoch B and then decays to flr^i + 62- Just for an illustration, let us choose 
61 = 10~^, Qr = 10~^ and S2 = 10~^. In this case, curvature perturbation is temporarily 
amplified to 10~^ and then decays to the observed value 10"^. 

It is worthwhile to mention here that the existence of more than a single curvaton is 
essential to make the curvature perturbations decay in time. One may consider a simpler 
situation where the inflaton generates large amplitude of the curvature perturbations by 
the standard mechanism and a curvaton that acquires little quantum fluctuations during 



final perturbation ~ fi^^i + ^; 
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inflation dominates the Universe at late time. In this case, attenuation of the curvature 
perturbation due to the curvaton dominance does not happen. Instead, the large amplitude 
of the curvature perturbations of the inflaton origin is taken over by the curvaton because 
the epoch when the curvaton starts its oscillations is modulated. We need at least two 
fields other than the inflaton to have a sensible model in which the temporal enhancement 
happens. 

In the following sections, we will give more quantitative discussions of this scenario 
and some interesting consequences. 



3 Decays of two curvatons and generation of pertur- 
bations 

3.1 Evolution of the background quantities 

In this subsection, we see how the background quantities evolve in time. The periods we 
focus on are from epochs A to D whose definitions are given in secj2j 

The background spacetime is the spatially flat Friedmann-Lemaitre-Robertson- Walker 
(FLRW) universe whose metric is given by 

ds"^ = a^ij]) {-dif + 6ijdx'dx^) = -dt^ + a'^{t)Sijdx'dx^ . (2) 

Here r] is the conformal time. As we mentioned in the last subsection, we can treat 
the dynamics of the curvatons as collections of non-relativistic particles. Therefore, the 
energy-momentum tensor for each curvaton can be written as 

= pzu'^zUz, (3) 
where Z runs 1 and 2 and is the four- velocity of Z-field with a normalization condition 

While the total energy-momentum tensor obeys the conservation law, each energy- 
momentum tensor is no longer conserved because curvaton fields decay into radiation 

^,Ti;^)u = QiA)., (4) 

where A now represents radiation, 1 or 2. The r.h.s. represents the transfer of the energy 
and momentum of the fluid. 

At the background order, Q(a)u for each fluid is given by 

(5(1)0 = aVipi, (5) 
Q(2)o = ar2P2, (6) 
Q{r)o = -aripi - aT2p2- (7) 

Because of the isotropy of the background spacetime, all the spatial components are zero. 
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Therefore, the evolution equations for the background quantities are given by 

p'l + SHpi = -alipi, pi + 3ifpi = -ripi, (8) 

P2 + 'iUp2 = -aT2p2, P2 + 3ifp2 = -T2P2, (9) 

p'r + AT-ipr = aTipi + aT2P2, Pr + 4:Hpr = Tipi + T2P2, (10) 

•H' = ^(Pl + P2 + Pr)«', i/^ ^ ^(pi + P2 + p,), (11) 

where a prime and an overdot denote differentiation with respect to 77 and t, respectively, 
with T-L = a' /a and H = a/a. From the first two equations, we have 

where is an arbitrary time in the epoch A. Substituting these solutions into the third 
equation, we get 

The designing situation in this paper is that both two curvatons dominate the universe 
at different epochs. We give in FigJT]a typical evolution of VLi and ^2- We see that their 
evolution can be clearly devided into four epochs explained in the last section. During 
epoch A, both of them grows as oc a. During epoch B, f2i ~ 1 and f22 ^ 1 stays constant. 
During epoch C, ^2 starts to grow again like oc a. During epoch D, ^2 — 1 and this epoch 
terminates by the decay of the (T2-field. 



3.2 Linear order perturbation equations 

Before discussing perturbation behaviors and their consequences, we give basic evolution 
equations for the scalar perturbations. 

We write the perturbed metric of the scalar type in the conformal Newtonian gauge: 

ds^ = a^{T]){ - (1 + 2i))drf + (1 - 2(t))5ijdx'dx^]. (15) 

For the matter perturbations, we intoduce density contrast and the velocity perturbation 

by 

OA = , UAi = av^. (16) 

Pa 

Note that Uq = — a(l + ip) is completely determined by the metric perturbation by using 
the normalization condition. 
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Figure 1: The evolution of Qi and Q2- We chose f2i,ini/fi2,ini = 10 and r2/ri = 10 
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At the hnear order in the perturbation, Q(^a)u for each fluid is given by 
SQ{i)Q = aTiSpi + aTipiip, 

SQ(r)o = —aTiSpi — aTipiip — aT2Sp2 — aT2P2ip, 

5Q{i)i = -aVipivl, 

SQ{2)i = -aT2p2v\, 

SQ{r)i = aVipivl + aT2p2V2- 



(17) 
(18) 
(19) 
(20) 
(21) 
(22) 



The second terms in f lT7|) and f|T8|) represent an effect due to the modulation of time. Q(r)^J. 
for the radiation is obtained from the conservation law for the total energy momentum. 
Using these equations, we flnd that the continuity equation and the Euler equation for 
each fluid are given by 



S[ — 30' + kvi = —aTiifj, 
62 — 30' + kv2 = —ar2ip, 

5; - 40' + hvr = (5i -Sr + ^p) + aT2— {62 - 5, + z^) 

3 Pr Pr 



(23) 
(24) 

(25) 



and 



v[ + Uvi -ktp = 0, 
f 2 + 'Hv2 — kip = 0, 

V„ - -Or - kip = all -Vi -Vr]+aL 2 -7V2 - Vr 

4 Pr \A J Pr \A 



(26) 
(27) 

(28) 
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respectively, where k is the comoving wavenumber. From the perturbed Einstein equations, 



we can derive the following evolution equations: 

+ 3n(j)' + 3n^ip = -4:7iG{pi6i + P2S2 + Pr5r)a'^, (29) 

kinip + (f)') = 4nGa^ (^piVi + ^2^2 + ^PrVr^ , (30) 

0-V^ = O. (31) 



3.3 Perturbation evolution on super-horizon scales 

Since understanding the evolution of the curvature perturbation on super-horizon scales 
is important for our purposes, let us next consider the perturbation behavior on super- 
horizon scales, postponing the analysis on sub-horizon scales of perturbations. The evolu- 
tion equations on super-horizon scales can be obtained by setting k = in the perturbation 
equations we gave in the last subsection. 

In FiglJl we show a typical evolution of (p which is obtained by solving numerically 
the perturbation equations with k = 0. We started the calculation at a time when the 
expansion rate H is much larger than Fi, with the initial conditions in the epoch A given 



by the analytic approximation: 

3 

Sliv) = ^l,mi - o ('^l,mi^l('7) + hhM2iv)) H > (32) 

o 

3 

^2(^7) = hini - o ('^l,ini^l('7) + ^2,ini^2(^)) H , (33) 

o 

1 

^riv) = -2 (^l,ini^^l(^) + ^2,inif^2(^)) + " " " , (34) 
<Piv) = -I (5l,ini^^l(r/) + 52,inif^2(^)) + " " " , (35) 



where ■ ■ ■ represent terms that are suppressed at early time with higher powers of 
and/or Since we assume that inflaton fluctuation contributes little to the curvature 
perturbation, we have imposed the condition that both 6r and vanish at the outset. The 
parameters used in Figj2] are such that 5i^ini = 10^^ and 52,im = 10"^. The parameters for 
the background are the same as the ones used in FigO 

From FiglH we can clearly see that is enhanced during epoch B. The order of (p at this 
time is roughly (5i^ini (we will provide exact analytic expression later.). This enhancement 
still persists during epoch C until the (T2-field dominates the universe. As the cr2-field 
dominates, (j) decays to f2r,2<5i,ini + <52,ini, where ^2^,2 ^ 1 is evaluated at the time 
when the (T2-field decays. Therefore, the enhanced value of the curvature perturbation is 
determined by Si^i^i and the final value is determined by either f2r-,2<^i,ini or (52,im, whichever 
is greater. 

We can analytically derive time evolution of curvature perturbations in the super- 
horizon regime solving perturbation equations for gauge-invariant variable, (a, which is 
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Figure 2: The evolution of cj). Since is negative, we plot —(f). We chose f^i,ini/^2,ini 

10, = 10-2, <52,ini = lo-^ r2/ri = 10-10. 



defined by = —0 — and represents curvature perturbation on the uniform en- 

ergy density surface of the component A. In the presence of interactions their evolution 
equations in the super-horizon {k — y 0) limit read 



C2 
Cr 



3 

Tipi 



r2C2 



(36) 
(37) 



[-(4 + 37i)Cr + (3 + 37i)Ci] + -f^ [-(4 + 372)Cr + (3 + 372)C2] , (38) 



wehre 'Ja = ^aH'^H). Using the background solutions for pi(t), P2it), and Pr(t), namely 



(fT2|) . (fT3|l . and (fT^ . these equations are solved as 

-^ri (t-t.) 





= Cue 


C2{t) 




Cr{t) 





X exp 



L + 7i)Ci.e-^^^(^'-') + -^^il + 72)C2*e-^^^(*'-*-) 
"^(4 + 370 + 4^(4 + 37^ 



dt" )■ dt' 



(39) 
(40) 



(41) 



In terms of these solutions the curvature perturbation in the uniform total density surface, 
C, is expressed as 

^ = _0 _ H— = ^^1^1 + ^^2(2 + 4pr-Cr + 37lPl(Cl - Cr) + 372P2(C2 " Cr) .^2) 

p 3pi + 3p2 + 4pr 

The key quantity to understand the time evolution of C is Cr- In the regimes A and B, only 
the terms involving pi is important in the both integrands in (BTj) because by assumption 
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P2 is much smaller than pi then. Furthermore, in the regime B when pi ^ the last 
exponential factor in the right-hand-side of ( 14T]) takes an appreciable value of order of 
unity only for a short time interval 

^t" = t-t'^ (43) 

ripi(4 + 37i) 

which also limits the range of t' integral. As a result we find 

We can estimate C,r at the end of the regime B by the contribution at the epoch when the 
ratio pi/ pr is the largest. This is just before cxi's decay when 71 was still negligible with 
Ciif) ~ Ci*- We therefore find Cr ~ 3Ci*/4 at the end of the regime B. 

The behavior of Cr in the regime D can also be understood similarly replacing the suffix 
1 by 2. When p2 dominates over p^, we find C^rif) ~ 3Ci*/4 + 3C2*/4, so that the total 
curvature perturbation reads 

^^^^ _ 3p2(t)C2* + 3p,.(t)(Cl. + C2) 

^'^ 3p2(t)+4p,(t) ' ^ ' 

which clearly shows that as the fraction p2{t)/ prit) increases the large curvature pertur- 
bation due to Ci* is regulated to a smaller value to reach the final value 



P2 



C,2* H 

d/ P2 



Here ^ 

P2 



/ N 2/3 n4/3 

I r7 ) I ^h^ ) denotes the ratio of radiation energy from pi to p2 upon 
decay of (J2. Note that Ci^* = 5i,*/3 and (2,* = 52,*/3 hold because is negligible and 
Fi, F2 <C at the outset. 



d 



3.4 Nonlinear super-horizon perturbations 

Eq. (l46l) clarifies how the final curvature perturbation after its temporal enhancement 
is related to the initial amplitudes of the two curvatons. Although an essential point is 
completely manifested by Eq. ( H6l) , the approximations we have made to derive it neglects 
a (9(1) constant factor in front of the second term in Eq. fH6|) . 

In this subsection, by using 6N formalism [391 - H3] . we derive the exact analytic expres- 
sions for the enhanced and the final amplitudes of the curvature perturbation not only to 
linear order in perturbation but also to any higher order. According to this formalism, 
the curvature perturbation on the uniform total energy density hypersurface at a point x 
is given by the perturbation of the e-folding number: 

C(?7,x) = A^(cri(?7, x), (72(77, x)) — spatial average, 

(ir/' 7/(77', 0"! (77', x), (72(77', x)) — spatial average, (47) 
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where 77* is an arbitrary earlier time than 77 and the hypersurface at 77* should be the 
flat slicing. The evolutions of c'"i,cr2 and "H are determined by solving the background 
equations fl8l)-f|TIl). 

In the appendix, we show that the number of e- folds from 77=,, at the epoch A to 7// well 
after the a2-field decay is given by 

1 r / u \ ^/"^ / u \ '^/'^ 



iV(ai(7//,f),a2(7//,x)) =-log| (1 + er) (^^"Y^J n^[^{r]*,x) + 1]2,*(^*.^) 

1 /3 

+ l,og|(£) a.Cr} + ilog^ (48) 

where er ~ 1.183 is a numerical value and p*/p/ is the total energy density at ?7*/?7/. cr 
is another numerical value. Eq. ( HHj) is obtained without using the sudden decay approx- 
imation which is frequently adopted in the literature. In the limiting case where both 
curvatons dominate the universe at different epochs separated far enough, (148!) is exact. 
Therefore, by combining ( HHl) with ( H71) . we can derive the exact expression of ( which is 
correct to any order in the perturbation. We can relate fii,*(7/^,, and Q2,*{,V*^^) with the 
density contrast for each field as 

fii,*('7*,^) = ^1,* [1 + 5i('7*,^)] , ri2,*('7*,^) = ^2,* [1 + ^2('7*,^)] • (49) 

If we choose 7/* to be deep in the epoch A, each density contrast can be approximated 
with the initial density contrast (see Eqs. (132 p and fl33|) ): 

^^l,*('7*, ^) = ^1,* [1 + 5l,ini(^)] , ^2,*{V*, ^) = ^2,* [1 + '^2,ini(^)] • (50) 

The initial density contrast is determined by quantum fluctuations generated when the 
mode crossed the Hubble length during inflation. 

Using these equations, we find that the full non-linear expression of ( is given by 

avf,x) =ilog |(1 + er)fif (7/.)[l + 6^Mx)f' + (^)'^'f^f (^*)[1 + S^Mx)^ 

-hoglil + er)n'/\v*)+(^Y'^Tiv*)]- (51) 



At linear order in the density contrasts, this equation reduces to 

«3(l + er)(^^J [j^J 5l,ini(x) + -52,ini(f). (52) 
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To get the second equation, we have used an inequahty <^ fig^^^ which is 

equivalent to the condition that the cr2-field dominates the Universe eventually. Therefore, 
it is useful to introduce a parameter s defined by 

s^il + eJ^Y'(^Y\ (53) 



rj 

which roughly represents the fraction of the radiation from ai-field decay at the time of 
(j2-field decay. With s, ( becomes 

C(?7/, X) ^ ^S(5i,ini(x) + ^52,ini(^), (54) 



which is a more regorous expression of (1461) . We see that the transfer coefficient of 6i is 
proportional to s. Apart from the numerical factors, Eq. (15^ is nothing more than the 
rough estimation ([1]) and also agrees with the result of |21] . 

We can also derive the maximum magnitude of the enhanced ( during the epochs B and 
C. Since the cr2-field is subdominant during those epochs, Cmax is completely sourced by the 
(Ti-field perturbation. This means that Cmax is equal to the final curvature perturbation in 
the single curvaton model in which the curvaton dominates the universe before its decay. 
Therefore, Cmax can be obtained by the second term in ( 15^ with 62 replaced by 61. 

Cmax(^) = ^5l,ini(^)- (55) 

We can convert C(^/) Cmax into the corresponding 0(?7/) and 0max- In the linear 
perturbation theory, it is well known that is related to C by [5] 

2 

= — -C, (for the radiation dominated universe) , (56) 
3 

(j) = — (, (for the matter dominated universe) . (57) 
5 

Therefore, 0(///) becomes 

2 2 

HVf) = -gS5l,ini - g^2,mi- (58) 

On the other hand, 0max at the epoch B/C is given by 



--^lini) (in the epoch B), (59) 
5 ' 

2 

-g^i.ini, (in the epoch C). (60) 
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4 Implications 



We have shown that the curvature perturbation can be temporarily enhanced in two 
curvaton models. At first glance, it seems that such an enhancement has nothing to do 
with observations since what we observe is the curvature perturbation at or after the time 
of last scattering epoch when it has already settled down to the observed value 10~^. This 
naive guess is true for the power spectrum. However, the trace of the enhancement enters 
the game when we consider the black hole formation, higher order correlation functions of 
the curvature perturbation (non-Gaussianity), and the gravitational waves generated by 
scalar-scalar mode couplings. 

4.1 Non-Gaussianity 

Eq. (15T|) is the fully non-linear expression of the curvature perturbation. By using this 
equation, we can calculate correlation functions of any order. In this paper, we calculate 
the three and four-point functions (bi- and tri-spectra) which are now becoming important 
observables to extract information of the early universe. 

The curvature perturbation given by Eq. ( ISTl) is the so-called local type for which the 
curvature perturbation depends on the source fields at the same point [S] 

C(x) = NaSaa{x) + l;NabSaa{x)6ab{x) + ^NabcS(Ta{x)Sab{x)6o-c{x) H , (61) 

z o 

where 6aa{x) is the Gaussian field fluctuation at a point x at some initial time and Na = 
dN/daa. 

The power spectrum, P^, bispectrum, B^, and trispectrum, T^, of the curvature per- 
turbation are defined by 

(Cg^Cg,) = (2vr)^Pc(^i)^(^i + ^2), (62) 
iCkSk^CkJ = i^^fBdh, /C2, h)5{h + h + 4), (63) 

and 

iCkSk,Q,Ck,) = (2vr)'Tc(fci, k2, h, h)5{ki + h + h + h), (64) 

respectively. For the case of local type curvature perturbation, P^ and can be written 

as 

Bcih, k2, h) = ^/nl (Pc(^i)n(^2) + Pc(^2)Pc(^3) + P(,{h)P^{h)) , (65) 

T^(fci, /C2, ^3, ki) = Tnl (P^(A;i3)P^(A;3)P^(A;4) + 11 perms.) 
54 

+7^9NL {Pdk2)Pdh)Pdh) + 3 perms.) , (66) 
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with kis = \ki + k^l. Here the constant parameters /nl,tnl and gjssh are the so-called 
non-linearity parameters and given by 



6 , NMN'" 



and 



54 iV^fe^iV^iV^iV^ , , 

^^^^ = (iV.iV^)^ • ^''^ 

Let us define a new parameter r by 

r = SO-i/o-2 ~ 52,im/(s5l,mi), (70) 

which represents the contribution of 52,ini to C compared to that of 5i,ini- If r = 0, then 
( is solely sourced by the fiuctuations in ai. If r ^ 1, then ( is mostly sourced by the 
fiuctuations in a2- With this parameter, the non-linearity parameters for two curvaton 
case are given by 

_ (15r^ + 80)r^. - 25 _ 25 

12(r^ + lf. -12(r2 + lf. + ^^'^' ^^'^ 
(9r^ + 112r2 + 64)r2s2 - SOr^s + 25 25 / 1 \ 

T"NL = 5 = 5 1" C ~ ! ( '2) 

4(^2 + 1)3^2 4(r2 + l)=^s2 

_ (225r^ + 3300r2)r2s2 - 1500r2s + 125 _ 125 O ( (73) 

~ 108(r2 + 1)^52 - 108(r2 + 1)^52 ^ UJ ' ^ 



Since s enters /nl in the denominator, s <^ 1, which is satisfied in the situation we are 
interested in, is a necessary condition to have large /nl- This condition can be qualitatively 
understood by expanding ( 154|) to second order in 6ai: 

3 




C=oM2-r^+{ — ) 1=C. + 4^C (74) 



where we have set ^2 = and Q = is the Gaussian part of (. We see that the 



second order coefficient which is, apart from the numerical factor, nothing more than /] 



NL 



contains an enhancement factor 1/s This mechanism to get large /nl is exactly the 
same as the single curvaton case in which /nl is inversely proportional to a fraction of the 
curvaton energy density at the time when it decays into the radiation [16]. 



"^^Strictly speaking, at second order perturbation, terms should also appear in (|74p . However, such 
terms only yield /nl of 0(1) and can be safely neglected when s <C 1. 
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Since /nl is bounded to be |/nl| ^ 100 from the observations [T], s cannot be smaller 
than 10~^ if ^2 = 0. Correspondingly, the observationally allowed maximum curvature 
perturbation when it is enhanced is at most 10~^/s ^ 10~^. However, things change when 
62 is also allowed to take non-zero amplitude. From ( I7T1) . we find that if fluctuations in 
(72 contribute more to the final curvature perturbation than the ai-field fluctuations, i. e. 
r > 1, then /nl is suppressed by a factor compared to a case with r = 0. Therefore, s 
smaller than 10~^ can satisfy the bound |/nl| ^ 100 if r is suitably chosen. In particular, 
if r > ^ 1, then /nl becomes as small as (9(1). At this level of /nl, non-linear 

evolutionary effects become important and it will not be easy to extract primordial /nl 
from observations. Interestingly, in this case, we find from ( |72l) and ( 1731) that tnl, S'nl ^ 
0{s~^^^) ^ 1. Therefore, strong non-Gaussianity appears in the trispectrum but not in 



the bispectrum In such a case, the trispectrum would be useful to search for non- 
Gaussianity. 

From ( 172]) and ( 175]) . we can also derive a unique relation between tnl and (^nl as 

Tnl 27 



fi'NL 5 v/nlJ ' ^ ^ 

which will be useful to observationally discriminate the two curvaton model from the other 
models that generate large non-Gaussianity. Consistency relations between /nl and qtsil, 
/nl and tnl or among the three parameters have been obtained for various models. As 
far as we know, ( 173]) is the first example that gives the unique relation between tnl and 
(j'NL without using /nl (for consistency relations for other models, see [17]). 



4.2 Stochastic gravitational waves 

It is well known that second order tensor perturbations are induced by the first order 
scalar perturbations by the mode- mode couplings. This means, ficw, the energy density 
of GWs per unit logarithmic interval of frequency, is proportional to quartic of the scalar 
perturbations, motivating us to consider the generation of GWs in the two-curvaton model. 

A basic picture of the production of GWs in the two curvaton model is that the scalar 
perturbations are temporarily enhanced equally on all the superhorizon scales and large 
magnitude of ^gw with a frequency equal to the Hubble parameter is produced at each 
time while the enhanced modes are re-entering the horizon. Therefore, we expect that the 
resulting Qqw will have a broad peak of an interval of frequencies whose corresponding 
modes re-enter the horizon during the curvature perturbation is being enhanced. To 
get this kind of result, we have to solve the equation of motion for the second order 

"^■^It is generally true that the trispectrum becomes relativily stronger than the bispectrum if more than 
one field contribute to the curvature perturbations. The local-type single field model yields a relation 
■''NL = H/nl- SS], it was shown that an inequality tnl > H/nl holds for any local- type multi-field 
model. In |49| . a scenario that both inflaton and curvaton contribute to the curvature perturbations was 
considered. It was shown that tnl is enhanced by a factor 1 + compared to where r represents 

the contribution of the inflaton fluctuations to the curvature perturbations. 
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tensor perturbations as well as the perturbation equations for the scalar modes and the 
background equations. Although all the numerical calculations to arrive at the final result 
are straightforward, it takes a long computation time due to multi-integrations. To avoid 
this, we stick ourselves to cases where ai-field immediately decays soon after it dominates 
the total energy density. In other words, we assume that the epoch B terminates in a 
moment and the effects of the enhancement of the curvature perturbation shows up only 
at the epoch C which is radiation dominated. In such cases, instead of numerically solving 
the perturbation equations for the scalar modes, we can assume that the scalar modes 
which are to re-enter the horizon during the epoch C are already enhanced since the epoch 
A and can use the analytic transfer function for the radiation dominated universe to evolve 
the scalar modes. By these assumptions, we do not need to numerically integrate the 
background equations and the linearized equations for the scalar modes, which drastically 
diminishes the task of the numerical computations. What is then left is to solve the 
evolution equations for the second order tensor perturbations sourced by the first order 
scalar perturbations. 

In what follows, we first briefiy review the general formalism to calculate ficw induced 
by the scalar-scalar couplings [22] - [2l] and then provide our results in two curvaton case. 
We basically follow the notations of [50]. The metric that is relevant to our purpose is 
given by the metric ( fT5l) plus the second order tensor perturbations: 

ds^ = a'^{r]) [ - (1 + 2tp)dri'^ + {(1 - 20)5^^ + hij}dx'dx^] , (76) 

where hij are the second order tensor perturbations which satisfy the transverse-traceless 
conditions: 

hijj = S'^h,j = 0. (77) 

In the following, we set if) = (j). Since the gravitational waves are transverse waves, we can 
Fourier-decompose hij as 

h^Av, ^) = / '''''' + -^^Ak)hM) , (78) 

where eij{k) and eij{k) are the polarization tensors orthogonal to k. Introducing the unit 
vectors ei{k) and ei{k) orthogonal to k, they are given by 

eij{k) = (^ei{k)ej{k) - ei{k)ej{k)^ , eij{k) = {ei{k)ej{k) + ei{k)ej{k)^ . (79) 

Then the evolution equation for hj: is given by 



hi + 2'Hhl + eh^ = S,;, (80) 



where the source term is given by 



^ r d^q ( 2 / 1 / / \ 

Si = '2e''{k) J q,qm \^<Pcfl>j:_^+ ^4>,fi^j:_^+ J^^'.^'k-g) ■ 
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In the radiation dominated universe, the solution of Eq. ( IHOj) is given by 

MV)= r dr^' guM)a{r,')Su{v'). (82) 



where = a/i^ and Qkiv^v') is the retarded Green's function: 

9k{v,v') = '^s^T^Kv -V')- (83) 

Then the two-point function of /ig becomes 
I rv rv 

(^fci^fc^) = y J dr]2 gkAv,Vi)9k2iv,V2)a{r]i)a{r]2){Sf:^{r]i)Sj^^{r]2)). (84) 

We need to evaluate (5'^^(?7i)5'^^(?72))- To this end, we need time evolution of </)^. Its 
evolution in the radiation dominated universe is given by 

0,,(r/) = Dkiv)Mk), Dkiv) = 77^ f ^ sin ( ^] - cos f ^ ) 1 . (85) 



Using this, the two-point function of the source can be written as 

^qiqjqmqnf{\ki - q\,q,vi)fi\ki - q\,q,V2), (86) 

where 

(00(^1)00(^2)) = {27rfP^{h)6{h + h), (87) 
f{h,k2,v) = 3DkMDkM+v{2DkM+vD'kM) D'kM- (88) 
We define the power spectrum of the tensor perturbation by 

= i2nfP,iv,k,)6ih + k2) (89) 

When the source term can be neglected, the energy density of the tensor perturbation can 
be written in terms of their power spectrum well inside the horizon, 

Uk^' 



now{v,k) = -i-] Vhiv.k), (90) 



where Vh = k'^Ph/i^ix) is the dimensionless power spectrum. 
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Figure 3: f^Gwl^;''?) for various A. The time is chosen to be 77 = 10^ /k. 

Given the general formahsm to calculate ficw? we are ready to apply it to the two 
curvaton model. Due to the reasoning we mentioned at the beginning of this subsection, 
we use top-hat type function with a width A for V,f,{k): 

V,{k) = [^ |ln(fc/MI<A, (91) 

I otherwise, 

where = k^P^/{2 tt) , kp is the peak wavenumber and A?' corresponds to the total power 
of the spectrum. Then VLqyj can be written as 



^Gw{v,k) = ^ I dq I d5 r\k,ri,q,5), (92) 

9min ^ —^n 



6 

where 

gmin = max{A;/2,e^'^}, (5max = minjl, ^ sinh A [>. (93) 

The function J-'(/c, r/, g, 5) is defined by 

k (V-A;2)(l-52) / \ 

where 

I{k,ki,k2,r]) = k dr]i a{r]i)gk{r],r]i)f{ki,k2,T]i). (95) 
Jo 

In Fig. [3], we show numerically calculated A~'^QGw{k,f])/{sQrad) for various values of 



A #1 A suppression factor s is due to the decay of ^gw during the epoch D where the 



^^The results presented here look somewhat different from the ones given in [5^ which also obtained 
f2GW for the top-hat type power spectrum, although the qualitative features are the same. This is due to 
there were a couple of mistakes in the manipulation in '50' . 
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Figure 4: Left panel shows contours of ficw.piatcau and the right one shows contours of 
/max and /min- The other free parameters A and Fi are fixed to be = 5 x 10~^ and 
Fi = 100 GeV. 



universe expands hke matter dominated regime and the other one Qrad — 8 x 10"^ comes 
from the matter dominated universe after the matter-radiation equahty. We see that for 
A < 0.5, has a strong peak at k = kp and decays as k~^ for k <^ kp. In all the 

cases, l^Gw has a sharp drop at /c/Zcp ~ e^. This is because the momentum conservation 
prohibits a generation of the tensor mode whose wavenumber is greater than 2kpe^. For 
A > 0.5, l^Gw has plateau between a range < k/kp < and decays as k~^ for 
k < e~^kp. The magnitude of the plateau is well fitted with ^ 2.0 x A~^. Therefore, the 
magnitude of the plateau at present is given by 

nGW,plateau - ^^^-ad - 2 X 10"^^ ^ (^) (96) 

The current frequency corresponding to the upper limit of the plateau is 

1/2 



^ V^^V4 ^ 10 Hz r^V'Y^^V (97) 
^'"^^ - 2vr^ea " VIO-V I 100 GeV ' ' 



where ^eq — 3200 and H^q are the redshift and the Hubble parameter at the time of the 
matter-radiation equality respectively. The lowest frequency is 



/min j ■f^^^' ('^^) 

It is interesting to derive an upper bound on f^GW.piateau • From the expression of /nl 
given by Eq. (jTTl) and the WMAP normalization [1], 

V^{Vf) = {l + r')s'^^W-^ (99) 
18 
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Figure 5: This panel shows contours of /nl and tnl- The other free parameters A and Fi 
are fixed to be ^ 5 x lO'^ and Fi = 100 GeV. 

we can write s and r in terms of /nl and A^/A. By using these relations, we find that 
^GW,piateau IS bouuded from above as 



where in the last inequality, we used A > 0.5, /nl ^ 50 and A < 0.05. As we will see in the 
next subsection, larger value than A = 0.05 over-produces black holes, which is excluded 
by the observations. Therefore, f^Qw.piatoau = 10~^^ is the possible maximal amplitude that 
can be achieved in principle. This amplitude is much smaller than the one given in |50] . 
This is mainly because ^cw.piatcau in the two curvaton model considered here is multiplied 
by the suppression factor s (see Eq. (196 p ) due to the existence of matter dominance epoch 
(epoch D). On the other hand, [SU] does not consider such an epoch and hence there is no 
additional suppression on Qqw- As can be seen from Fig. |6l ficw.piateau = 10~^^ is slightly 
lower than the sensitivity of LISA |1] and DECIGO/BBO [52l[M]. But it is much higher 
than the sensitivities achieved by ultimate-DECIGO and space-based AGIS [51] . 

As an example, we show contour plots of /^ax, /mm and r^cw.piateau in Fig. Hand /nl 
and Tnl in Fig. O as functions of r defined by Eq. fl70]l and Q2,*/^i,*- The other free 
parameters A and Fi are fixed to be = 5 x 10~^ and Fi = 100 GeV. A is determined 
by a relation 



and s is determined by the WMAP normalization fl99l) . 

We find that f^cw.piateau is O(10~^^) in the frequency band 10~^ Hz — 30 Hz. We also 
find that in this case, r < 1 is ruled out by the observational bound on /nl, |/nl| < 100. 
This is because s is typically 10"'^ in this case and we have /nl = O(IO^) for r < 1. As we 




(100) 




(101) 
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Figure 6: This panel show plots of ^gw for three cases: A = 1.0, 3.0, 5.0, assuming 
Fi = 100 GeV, = 5 X 10"'* and r = 3 with expected sensitivity of the future GW 
detectors such as LISA, DECIGO/BBO and AGIS. 



have seen in the previous subsection, tnl becomes very large as O(IO^) — (9(10"^). Hence 
the strong non-Gaussian signal appears in the trispectrum rather than in the bispectrum. 

In Fig. [6l we show plots of flew for three cases: A = 1.0, 3.0, 5.0, assuming Fi = 
100 GeV, = 5 X 10"^ and r = 3 with expected sensitivity of the future GW detectors 
such as LISA, DECIGO/BBO and AGIS. As mentioned before, we see that maximum of 
flaw is below the LISA and DECIGO/BBO sensitivities. On the other hand, ultimate- 
DECIGO and space-based AGIS will be able to probe GWs generated in two curvaton 
models if the model parameters are suitably chosen. If we lower Fi as small as 0.1 meV, 
then ^GW.piateau does not change but /max becomes O(10~^ Hz). Hence ^cw.piatcau enters 
the frequency region of the pulsar timing. However, as is clear from Fig. [HI f^cw.piatcau is 
far below the upper limit coming from the pulsar timing observations. 



4.3 Primordial black holes 

PBHs are formed if there exist density fluctuations of the order unity and when such 
modes re-enter the Hubble radius [25l[27]. The mass of the PBHs is roughly equal to the 
horizon mass at the time of formation, 

GMpBH ^ ^. (102) 
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To simplify the analysis and also to make it consistent with the study of the last sub- 
section, we again assume that the epoch B terminates in a moment. Then, PBHs would 
be efficiently formed during the epoch C and the mass range of the resultant PBHs is 
estimated as 



r,<GA/pB„<f^^^j. (103) 

Let us write the fraction of the energy density of PBHs of mass between (M, M + dM) 
at the time of formation as 

(104) 

dM ^ ' 

Then the total fraction of the energy density of PBHs is given by the integral of the above 
quantity, 

/3PBH = \dM (105) 

where the range of integration is given by Eq. f llOSp . The purpose of this subsection is 
to provide /3pbh for the two curvaton model and to discuss its cosmological implications. 
Since most of the necessary formulae and their detailed derivations are written in [50], we 
omit the intermediate calculations and provide only results. 

Assuming that the gravitational potential smoothed over the horizon size is Gaussian, 
dii^^YildM can be written as 

where Rm is the comoving horizon length, i.e. qRm = GM, with a being the scale factor 
at the time of horizon crossing, 0c is the threshold value of for black hole formation and 

is the variance of smoothed over the horizon size. Although 0c depends on the initial 
configuration of the perturbations [55l[56], we simply use 0c = 0.5 [57] . 

For the top-hat type of V^{k) which we have considered in this paper, a^j^^ becomes 

al =:^logf-^y (108) 

Having these, we are ready to do the integration (11051) . The result is 

^,,„../^Mexp(-iLV (109) 



It A "V 2.42, 

Fig. [7] shows fipBH^^ as a function of A. A is chosen to be A = 1.0. It is clear that 
/5pBH is very sensitive to change of A. Since various cosmological observations put different 
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Figure 7: This panel shows /3pbh as a function of A. A is chosen to be A = 1. 

constraints on the abundance of PBHs for different mass, using Ppbu to constrain A is not 
correct in a precise sense. However, because of the strong sensitivity of /3pbh to A, even 
the change of many order of magnitude of /3pbh corresponds to the change of a factor of 
a few in A. This fact vahdates the use of /Spbh to constrain ^ as a first approximation. 
According to [28], the constraint on /3pbh varies from 10"^'' to 10"^''. This can be converted 
to a constraint on ^ as ^ < 0.05. 

Since ctrj,^ is a monotonically decreasing function of M, the mass of PBHs that domi- 
nantly contribute to (11091) is the lower limit of Eq. (I103p . which is estimated as 



Combining this with Eq. ( 1971) gives a relation between BH mass and the corresponding 
frequency of GWs, 



5 Summary 

If more than one curvaton dominates the Universe at different epochs, the curvature 
perturbations can be temporarily enhanced to a value much larger than the observed one 
10~^. In this paper, we studied in detail the evolution of the curvature perturbations in two 
curvaton models. By solving the linearized perturbation equations both numerically and 
analytically, we confirmed that the curvature perturbations are indeed enhanced during 
the period from the time when the first decaying curvaton dominates the Universe until 
the second decaying curvaton dominates the Universe. The amplitude of the enhanced 




(110) 





22 



curvature perturbation is roughly equal to the density perturbation of the first decaying 
curvaton. We then provided an analytic expression of the full-order curvature perturbation 
which does not rely on the sudden decay approximation and is exact on super-horizon 
scales. At the linear order, we compared with the analytic results with the numerical ones 
and found they agree very well. 

The temporal enhancement of the curvature perturbations leaves its traces as the strong 
non-Gaussian perturbations, stochastic gravitational waves generated by the scalar-scalar 
mode couplings and the primordial black holes. By using the analytic formula for the 
full-order curvature perturbations, we gave the expressions for the so-called non-linearity 
parameters /nl, tnl and ^^nl- If both two curvatons contribute to the final curvature 
perturbations, then the strongest non-Gaussian signal comes from tnl and (^nl rather than 
from /nl- We also gave a consistency relation between tnl and (7nl without using /nl- 
If non-Gaussianity is detected in the future, this relation would be useful to discriminate 
this model from the others that also generate large non-Gaussian perturbations. 

We next studied the generation of GWs sourced by the enhanced curvature pertur- 
bations. The spectrum of Qqw has a plateau corresponding to the duration of the en- 
hancement. Because of the existence of the period where the secondly decaying curvaton 
dominates the Universe and the Universe expands like the matter dominated universe, 
^GW today is accompanied by a suppression factor which represents the fraction of the 
radiation energy density coming from the first decaying curvaton at the time when the 
secondly decaying curvaton decays. Due to this suppression factor, the possible maxi- 
mal amplitude of ^gw at the plateau is at most 10"^^, which is below the LISA and 
DECIGO/BBO sensitivities, but above the ultimate-DECIGO and space-based AGIS sen- 
sitivities. Actually, if the decay rate of the first decaying curvaton is around 100 GeV, 
then the frequency interval of the plateau can be 10~^ Hz ~ 10 Hz. Such a case can be a 
target of the uhimate-DECIGO and AGIS. 

We finally calculated the abundance of PBHs which are formed by the gravitational 
collapse of the enhanced curvature perturbations. We then provided an upper bound on 
the amplitude of the enhanced curvature perturbation by using the observational upper 
bounds on the abundance of PBHs. 
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A Derivation of Eq. (1481) 



In this appendix, we derive Eq. fl48p . In the following, we frequently use epoch A, epoch 
B, epoch C and epoch D which are defined in SecHJ For convenience, we use the cosmic 



23 



time t as a time variable. The basic equations we use are the background equations ([8]), 
(Q for the curvatons and (ITOl) for the radiation. In the spirit of the 6N formahsm, the 
solutions of these equations, (fT2l) . (fT3l) . and ( fT4l) depend on the space coordinate x and 
evolve independently of other points, although we do not show their dependence explicitly. 
Also, t should be interpreted as the local proper time at x hereafter. 

If we choose t such that ^ 1, z.e., an epoch well after the (j2-field decay, then p^ad 
can be safely replaced by the total energy density p{t). Then from f|T^ we find 

a\t)p{t) = dt' {V^a{t')alpi^,e-^^'' + T2a{t')alp2,,e-^^''^ . (112) 

Also since the integrand is exponentially suppressed for t ^ l/r2, we can push the upper 
limit of integration to infinity. Then, the e-folding number from the initial time to the 
time when the total energy density becomes Pf is given by 

= ^ log (rifii,,Fi + r2fi2,*F2) + \ log (^^^ , (113) 

where Fi and F2 are defined by 

/■oo 

Fi= I dta{t)e-^^\ (114) 

fOO 

dta{t)e-^^K (115) 





We see that the e-folding number depends on the field values of a and h through Fi and F2 
as well as and 1^2,*. Now the problem is reduced to deriving the analytic expression 
of Fi and F2. 

Let us first evaluate Fi. To evaluate it, we need to know the evolution of the scale factor. 
Since the integrand of Fi has an exponential cutoff for t^l /Fi, the dominant contribution 
to the integral comes from a time interval where ai-field dominates the universe (epoch B). 
During that epoch, we can neglect the (T2-field because its energy density is tiny compared 
to the total one. Then, the background equations in terms of VLi and H are given by 

^ = {(i-f!i)i/-r,}fi,, (116) 

f ^4(4-.,)//^. (117) 

We want to solve these equations from the time t = tb when the ai-field is dominating the 
universe, but still well before the cxi-field decay, to the time well after the cxi-field decay. 
To know the magnitude of the scale factor at t^, we need to connect the scale factor in 
that epoch with the initial scale factor a*. We can make the connection by solving the 
differential equation for a{t): 

H^ = Hl{[^\ (118) 
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We can exactly integrate this equation. The result is 



2 ^2 - 3^1,, + Jl+(—-l] n,^, (-2+ (— + 2] fii,, ) )■ = H,{t - Q. (119) 



Neglecting t*, we find that a(t) at tb is given by 

a{h) = a. (^^fii,,^ ^'{HJbf' + _^ + O (^1^ . (120) 

For large tb, the first term dominates the scale factor. The second term represents the 
contribution from the radiation being diluted faster than the cxi-field. As we will see, 
keeping only the first term gives the desired expression of Fi. 

Things become clearer if we introduce the new dimensionless variables by 

b{t) = -, (121) 

m ^ jH{t)- (122) 

Defining the dimensionless time s by s = Tit, the evolution equations for the new variables 
are given by 

(.3) 

f 4{l-(4-<}, (124) 

^={?^(1-SJ,)-1}S1,^ (125) 

as ^ 3s ' 

The initial conditions are 6(0) = -E(O) = f2i(0) = 1. The new set of differential equa- 
tions plus the initial conditions are free of model parameters, which means that we have 
extracted the dependence of Fx on the model parameters. Indeed, in terms of the new 
variables. Fx can be written as 

n \ 1/3 



F. = a.\-^^x,Hi^ (126) 

where cr is a purely numerical value defined by 

cr= ds s^'%{s)e-' . (127) 
We find numerically that cr ~ 0.830. 
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Let us next evaluate F2. For F2, the dominant contribution to the integral comes from 
a time interval when the (T2-field dominates the universe (epoch D). Since ai-field has 
completely decayed into the radiation by that epoch, we can only consider the radiation 
and the cr2-field in the following analysis. To make an argument similar to the case of Fi, 
first we have to connect the scale factor at the epoch C with a*. The scale factor at the 
epoch C is given by b{s) for s ^ 1. From fll23p and f ll24p . we find that the asymptotic 
form of b{s) for s ^ 1 is given by 

b{s) = drs-^/^ (1 + 0{s-^)) , (128) 

where dr is a numerical constant. We find numerically dr ~ 1.101. Substituting this into 
f ll2ip . the scale factor at a time t = tc itl the epoch C is given by 

a{Q = a, dr^fj'. (129) 

Using this result, the fl2 at this time is found to be 

Now it may be expected that the things can be exactly mapped to what we did for 
Fi by replacing the epoch A/B by the epoch C/D. To be more precise, simply replacing 
a^., f2i and Fi appearing in f ll26p by a(tc) given by fll29p . ^2 given by fllSOp . H{tc) 
and F2 would give the desired expression for F2. However, as it will come out, we need to 
include a term in the scale factor corresponding to the second term of fll20p to have the 
correct answer, which makes the calculations more complicated. 

If we denote by the time when the (T2-field is dominating the universe, but still well 
before the (j2-field decay, the scale factor at that time is given by 

aiU) = a{U) {^-n2{QH\tS^ "\f + ^ + ^ (^) ' (131) 

which is obtained by exactly the same argument as what we used to derive fll20p . After 
the example of Fi, let us introduce the new dimensionless variables by 



d{t) - "^1, ..vs.,. - (132) 



i{Q{p2{tc)H^{tc)y'"ty^ 

To include the second term of fll3ip . let us decompose the scale factor as 

a{t) = d{t) + Sa{t). (133) 

Here d{t){6a{t)) is a part of the scale factor which reduces to the first(second) term in 
(11311) at td- Correspondingly, we can define d{t), 6d{t), E{t), SE{t), ri2{t) and 6fl2{t)- 
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The evolution equations for d(t), E(t) and fl2(t) are given by fll23l) - fll25l) with b, E, Qi 
and s being replaced by d, E, Qi and u = r2t. The initial conditions for those variables 
are d{0) = E{0) = ^2(0) = 1. 

The evolution equations for 6d, 6E and 6Q2 at linear order are obtained by perturbing 
flT23D-ffT25ll: 

—5d = ^- ^d + —SE, 134 

au 6U 6U 

'^-5E = 1(1 - (4 - ^2)^\5E + -{- ^^5E + (135) 



du uy 3j n|^3 3^ 

= ((1 - 0^)^ - l| + n^l^-^^^SE - ^5^2). (136) 
du y 3m J [ 3m 3m J 

Behaviors of 5d, 6E and 6Q2 for m <^ 1 are obtained as follows. Combining 013ip and 
f ll32p yields the time evolution of 6d{s) for w <C 1 as 

p2/3 

Sd{u) = '-^^ . (137) 

(ffi2(tc)i/2(tc)) ^ f^2(tc)w2/3 

For 6E, by using a relation 6E = ^SH, we find that 

p2/3 

6E{u) = ^— . (138) 

(|fi2(tc)i/2(t,))^/'fi2(tc)«2/3 

For (5^2; since it is sourced by the residual radiation, we have 

oihiu) = — -— = — — ,^ , , = -7:5 . (139) 

Pi(t) a{u)Q2{tc) {\n2{t,)H\t,)y"'n2{t,)u^/^ 

To simplify the system further, let us introduce the new variables by 

1/3 



X ^ {hl2{QH\tS) (140) 
Y ^ {^^n2{QH\tS) '''^^u'"5E, (141) 
Z ^ {^n2{QH\tS) '''^^n'"5n2. (142) 
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Then the evolution equations for these variables become 



>^ + - |(4 - 1^2)}>^ + ^ (-(4 - n2)Y + EZ) , (144) 
Z + 1^(1 - O2) - ijz + ^ ((1 - n,)Y - EZ) . (145) 



The corresponding initial conditions are X{0) = —Y{0) = —Z{0) = 1. It is now clear 
that the new set of differential equations plus the initial conditions are free of model 
parameters, which means that we have extracted the dependence of F2 on the model 
parameters. Indeed, we find that F2 can be written as 

where /r is a constant defined by 

poo 

fr= rfMX(M)e-" ^ 1.18716. (147) 
Jo 

The second term in fll46p comes from 6a. Substituting f ll26p and ( I146p into (IllSp . we get 



4 t ' ' V Ti y '* V r2 y } 4 " {\4: J j 4 ° p/ 

(148) 

where er is a numerical constant defined by 

Numerically, we find that er ~ 1.183. Since the last two terms are merely constants, only 
the first term contributes to the curvature perturbation. 
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